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Abstract 



Let F[d; a, 6] be the ring of Ore polynomials over a field (or skew field) F, where a is a automor- 
phism of F and 8 is a a-derivation. Given a matrix A 6 F[d; o, S]" xn , we show how to compute 
the Hermite form H of A and a unimodular matrix U such that UA — H. The algorithm requires 
a polynomial number of operations in F in terms of n and the degrees (in d) of the entries in 
A. When F = k(z) for some field k, it also requires time polynomial in the degree in z, and if 
k = Q it requires time polynomial in the bit length of the coefficients as well. Explicit analyses 
are provided for the complexity, in particular for the important cases of differential and shift 
polynomials over Q(z). To accomplish our algorithm, we develop the Dieudonne determinant 
and quasideterminant theory for Ore polynomial rings to get explicit bounds on the degrees and 
sizes of entries in H and U. 



1. Introduction 

The Ore polynomials are a natural algebraic structure which capture difference, q- 
difference, differential, and other non-commutative polynomial rings. Linear algebra over 
these rings are important in solving the corresponding rings of differential and difference 
equations. The basic concepts of pseudo-linear algebra is introduced nicely in Bronstein 
and Pctkovsek (1996); see Ore (1931) for the seminal introduction. 

On the other hand, canonical forms of matrices over commutative principal ideal 
domains (such as Z or F[x], for a field F) have proven invaluable for both mathematical 
and computational purposes. One of the successes of computer algebra over the past 
three decades has been the development of fast algorithms for computing these canonical 
forms. These include triangular forms such as the Hermite form (Hermite, 1863), low 



Email addresses: mwgSuwaterloo . ca (Mark Giesbrecht), ms2kim@uwaterloo . ca (Myung Sub Kim). 



Preprint submitted to Elsevier 



19 September 2011 



degree forms like the Popov form (Popov, 1972), as well as the diagonal Smith form 
(Smith, 1861). 

Canonical forms of matrices over non-commutative domains, especially rings of differ- 
ential and difference operators, are also extremely useful. These have been examined at 
least since Dickson (1923), Wcddcrburn (1932), and Jacobson (1943). Recently they have 
found uses in control systems (Chyzak, Quadrat, and Robcrtz, 2005; Zerz, 2006; Halas, 
2008). Computations with multidimensional linear systems over Ore algebras are nicely 
developed in Chyzak, Quadrat, and Robcrtz (2007), and a excellent implementation of 
many fundamental algorithms is provided in the OreModules package of Maple. 

In this paper we consider canonical forms of matrices of Ore polynomials over a skew 
field F. Let a : F — > F be an automorphism of F and 5 : F — > F be a er-derivation. 
That is, for any a, b G F, 5 (a + b) = 5(a) + 5(b) and 5(ab) = a (a) 5(b) + 5(a)b. We then 
define F[<9; a, 5] as the set of usual polynomials in F[d] under the usual addition, but with 
multiplication defined by 

da = <r(a)d + 5(a) 

for any a € F. This is well-known to be a left (and right) principal ideal domain, with a 
straightforward euclidean algorithm (see Ore (1933)). 

Some important cases over the field of rational functions F = k(z) over a field k are as 
follows: 

(1) <t(z) — S(z) — z + 1 is the so-called shift automorphism of k(z), and <5 identically 
zero on k. Then k(z)[d;S, 0] is generally referred to as the ring of shift polynomials. 
With a slight abuse of notation we write k(z)[<9;6>] for this ring. 

(2) 5(z) = 1 and a(z) = z, so 5(h(z)) = h'(z) for any h £ k(z) with hi its usual derivative. 
Then k(z)[<9; <r, 5] is called the ring of differential polynomials. With a slight abuse of 
notation we write k(z)[<9;'] for this ring. 

A primary motivation in the definition of k(z) [d; '] is that there is a natural ac- 
tion on the space of infinitely diffcrentiable functions in z, namely the differential 
polynomial 

a m d m + a^d™- 1 + • • • + aid + a e k(z)[d;>] 
acts as the linear differential operator 

on a differentiable function f(z). See Bronstein and Petkovsek (1996). Solving and 
analyzing systems of such operators involves working with matrices over k(z) [(?;'], 
and invariants such as the differential analogues of the Smith, Popov and Hermite 
forms provide important structural information. 
A matrix H 6 F[<9; er, 5] n n of full left row rank is said to be in Hermite form if H 
is upper triangular, if every diagonal entry is monic, and every off-diagonal entry has 
degree less than the diagonal entry below it. For example, in the differential polynomials 
<Q(z)[d; '] as above: 



A 



( i + (z+2)d + d 2 2 + (2z+i)<9 i + (i+z)d ^ 

(2z+z 2 ) + zd (2+2z+2z 2 ) + 8 lz+z 2 

^(3+z) + (3+z)d + d 2 (8+4z) + (5+3z)d + 3 2 (7+8z) + (2+4z)<9 j 



£Q(z){d;'] 3x3 (1.1) 



2 



has Hcrmitc form 



H 



f(2+z)+d l+2z -2+z+2z 2 

(2+z) + d i+ 7 -f + \d 
\ -f + =M^+*l d + d 2J 



6QW19; 



. /I3x3 



Note that the Hermite form may have denominators in z. Also, while this example does 
not demonstrate it, it is common that the degrees in the Hermite form, in both z and d, 
are substantially larger than in the input. 

For any matrix A 6 F[<9; er, S] n n of full left row rank, there exists a unique unimodular 
matrix U G F[d;a, S] n ™ (i.e., a matrix whose inverse exists and is also in F[d;a, S] n n ) 
such that UA — H is in Hermite form. This form is canonical in the sense that if two 
matrices A, B £ F[<9;er, S] n " are such that A = PB for unimodular P e F[d; a, 5] n ™ 
then the Hermite form of A equals the Hermite form of B. The existence and uniqueness 
of this form over F[d; cr,S] n " is established much as it is over Z™ x ™, with the division with 
remainder and the euclidean algorithm replaced by right division with remainder and the 
right euclidean algorithm respectively. See, for example, Newman (1972), Theorem II. 2. 

In commutative domains such as Z and F[x] there have been enormous advances in 
the past two decades in computing Hermite, Smith and Popov forms. Polynomial-time 
algorithms for the Smith and Hermite forms over F[x] were developed by Kannan (1985), 
with important advances by Kaltofcn, Krishnamoorthy, and Saunders (1987), Villard 
(1995), Mulders and Storjohann (2003), and many others. One of the key features of this 
recent work in computing canonical forms has been a careful analysis of the complexity in 
terms of matrix size, entry degree, and coefficient swell. Clearly identifying and analyzing 
the cost in terms of all these parameters has led to a dramatic drop in both theoretical 
and practical complexity. 

Computing the classical Smith and Hermite forms of matrices over Ore domains has 
received less attention though canonical forms of differential polynomial matrices have 
applications in solving differential systems and control theory (see Halas (2008); Kotta, 
Leiback, and Halas (2008)). Abramov and Bronstein (2001) analyzes the number of re- 
duction steps necessary to compute a row-reduced form, while Bcckermann, Cheng, and 
Labahn (2006) analyze the complexity of row reduction in terms of matrix size, degree and 
the sizes of the coefficients of some shifts of the input matrix. Beckermann et al. (2006) 
demonstrates tight bounds on the degree and coefficient sizes of the output, which we will 
employ here. For the Popov form, Cheng (2003) gives an algorithm for matrices of shift 
polynomials. Cheng's approach involves order bases computation in order to eliminate 
lower order terms of Ore polynomial matrices. A main contribution of Cheng (2003) is 
to give an algorithm computing the left row rank and a row-reduced basis of the left 
nullspace of a matrix of Ore polynomials in a fraction-free way. This idea is extended in 
Davies, Cheng, and Labahn (2008) to compute Popov form of general Ore polynomial 
matrices. They reduce the problem of computing Popov form to a nullspace computa- 
tion. However, though Popov form is useful for rewriting high order terms with respect 
to low order terms, we want a different canonical form more suited to solving system of 
linear diophantine equations. Since the Hermite form is upper triangular, it meets this 
goal nicely, not to mention the fact that it is a "classical" canonical form. An implemen- 
tation of the basic (exponential-time) Hermite algorithm is provided by Culiancz (2005). 
In Giesbrecht and Kim (2009), we present a polynomial-time algorithm for the Hermite 
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form over Q(t)[d;']. While it relies on similar techniques as this current paper, the cost 
of the algorithm is considerably higher, and the coefficient bounds weaker, as well as not 
obviously working for all Ore polynomials. 

The related "two-sided" problem of computing the Jacobson (non-commutative Smith) 
canonical form has also been recently considered. Blinkov, Cid, Gcrdt, Plcsken, and 
Robertz (2003) implement the standard algorithm in Janet. Levandovskyy and Schindelar 
(2011) provide a very complete implementation, for the full Ore case over skew fields, of 
a Jacobson form algorithm using Grobner bases in Singular. Middckc (2008) has recently 
demonstrated for the Jacobson form of a matrix of differential polynomials, which requires 
time polynomial in the matrix size and degree (but the coefficient size is not analyzed). 

One of the primary difficulties in both developing efficient algorithms for matrices 
of Ore polynomials, and in their analysis, is the lack of a standard determinant, and 
the important bounds this provides on degrees in eliminations. In Section 2 we establish 
bounds on the degrees of entries in the inverse of a matrix over any non-commutative field 
with a reasonable degree function. We do this by introducing the quasideterminant of 
Gel'fand and Retakh (1991, 1992) and analyzing its interaction with the degree function. 
We also prove similar bounds on the degree of the Dieudonne determinant. In both cases, 
the bounds are essentially the same as for matrices over a commutative function field. 

In Section 3 we consider matrices over the Ore polynomials and bound the degrees of 
entries in the Hermite form and corresponding unimodular transformation matrices. We 
also bound the degrees of the Dieudonne determinants of these matrices. 

In Section 4 we give an algorithm that, given a matrix A £ F[d; a, S] nxn (of full left row 
rank), computes H and U such that UA = H, which requires a polynomial number of 
operations in n, deg A, and the size of the coefficients of the entries of A. In the common 
case where F = Q(z) these algorithms require time polynomial in n, deg d (A), deg z (A), 
and the bit-length of the rational polynomial coefficients. Our algorithms could no doubt 
be generalized to the non-square and non-full-rank cases, but we do not pursue this here. 
We also give explicit upper bounds on the degrees and coefficient sizes of U and H. 

2. Non-commutative determinants and degree bounds for linear equations 

One of the main difficulties in matrix computations in skew (non-commutative) fields, 
and a primary difference with the commutative case, is the lack of the usual determinant. 
In particular, the determinant allows us to bound the degrees of solutions to systems of 
equations, the size of the inverse or other decomposition, not to mention the degrees 
at intermediate steps of computations, through Hadamard-like formulas and Cramer's 
rules. 

The most common non-commutative determinant was defined by Dieudonne (1943), 
and is commonly called the Dieudonne determinant. It preserves some of the multiplica- 
tive properties of the usual commutative determinant, but is insufficient to establish the 
degree bounds we require (amongst other inadequacies). In Gel'fand and Retakh (1991, 
1992) introduced quasideterminants and a rich associated theory as a central tool in lin- 
ear algebra over non-commutative rings. Quasideterminants are more akin to the (inverse 
of the) entries of the classical adjoint of a matrix than a true determinant. We employ 
quasideterminants here to establish bounds on the degree of the entries in the inverse of 
a matrix, and on the Dieudonne determinant in this section, and on the Hermite form 
and its multiplier matrices in Section 3 
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We will establish bounds on degrees of quasideterminants and Dieudonne determinants 
for a general skew field K with a degree junction deg : K -> ZU {—00} satisfying the 
following. For a,b £ K: 

(i) If a ^ then deg a £ Z, and degO = —00; 

(ii) deg(a + b) < maxjdeg a, deg b}; 

(iii) deg(afe) = deg a + deg b; 

(iv) If a 7^ then deg(a~ 1 ) = — deg a. 

For example, if K = F(z) for some commutative field F and indeterminate t, for any 
a e K* we can define dega = degajv — dega^i, where aN,aD € F[:r] are such that 
a = ajv/a_D- 

2.1. Quasideterminants and degree bounds 

Following Gel'fand and Rctakh (1991, 1992), we define the quasideterminant as a 
collection of n 2 functions from K rex " KU {!}, where _L represents the function being 
undefined. Let A £ K nxrl and p, q € {1, . . . , n}. Assume A pq £ K is the (p, q) entry of K, 
and let A^ £ k("- 1 ) x ("- 1 ) be the matrix A with thepth row and qth column removed. 
Define the (p, <7)-quasidctcrminant of A as 

\A\ pq — A pq — ^ ApidA^lji) 1 Aj q , 

where the sum is taken over all summands where is defined. If all summands 

have l-A^^ljj undefined then \A\ pq is undefined (and has value _L). See Gel'fand and 
Retakh (1992). 

Fact 2.1 (Gel'fand and Retakh (1991), Theorem 1.6). Let A £ K nxn over a (possibly 
skew) field K. 

(1) The inverse matrix B = A^ 1 £ K" xn exists if and only if the following are true: 

(a) If the quasideterminant \A\u is defined then \A\ij 7^ 0, for all i,j £ {1, . . . , n}; 

(b ) For all p £ {1, . . . , n} there exists a q £ {1, . . . ,n\, such that the quasidetermi- 
nant \A\ pq is defined; 

(c) For all q £ {1, ... ,71} there exists a p £ {1, . . . , n} such that the quasidetermi- 
nant \A\ P q is defined; 

(2) If the inverse B exists, then for i,j £ {1, . . . , n} we have 

_\(\A\ ij y 1 if \A\ij is defined, 
1 if \A\ij is not defined. 

We now bound the size of the quasideterminants in terms of the size of the entries of 
A. Assume that K has a degree function as above. 

Theorem 2.2. Let A £ K nxn , such that either A {j = or < deg ,4^ < d for all 
i,j £ {1, . . . , n}. For all p,q £ {1, . . . , n} such that \A\ pq is defined we have — (n — l)d < 
\A\ pq < nd. 

Proof. We proceed by induction on n. 
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For n = 1, p = q = 1 and |^4|n = An, so clearly the property holds. Assume the 
statement is true for dimension n — 1. Then 



degL4| M = degU M - £ A^A^yy 1 A, 

where the sum is over all defined summands. Then using the inductive hypothesis we 
have 

deg \A\ pq < max <^ deg A, max { deg A pi - deg \A {pq) \ j{ + deg A jq \ \ 



and 



< 2d+ (n-2)d < nd, 
deg\A\ pq >-deg\A ( ^\ Jt > -(n-l)d. □ 



Corollary 2.3. Let A £ K nx ™ be unimodular, and B £ K nx ™ such that AB = I. Assume 
Aij = or < deg A v] < d for all i,j G {1, .. . , n}. T/ien degi? < (n - 

Proof. From Fact 2.1 we know that Bji = when is defined (and -Bji = 

otherwise). Thus deg Bji = — deg \A\ij < (n — l)d, and By = or degB^ > since A is 
unimodular. □ 



2.2. Dieudonne Determinants 

Let [K*,K*] be the commutator subgroup of the multiplicative group K* of K, the 
(normal) subgroup of K* generated by all pairs of elements of the form a~ 1 b~ 1 ab for 
a, b e K*. Thus K*/[K*, K*] is a commutative group. 

Let A £ K" x ™ be a matrix with a right inverse. The Bruhat Normal Form of A is 
a decomposition A = TDPV, where P £ K nx " is a permutation matrix inducing the 
permutation a : {1, . . . , n} — > {1, . . . , n}, and T,D,V £ K nxn are 



T = 



1 * 
1 



... *\ 



D = diag(ui, . . .,«„), 



V = 



1 o 



* i 



... o\ 



v* ••• 



V 



The Bruhat decomposition arises from gaussian elimination, much as the LUP decom- 
position does in the commutative case. We then define Set {A) — sign(er) ■ u\ ■ ■ ■ u n £ K 
(sometimes called the pre- determinant of A). Let 7r be the canonical projection from 
K* K/[K*, K*]. Then the Dieudonne determinant is defined as Pet (A) = tt(Sst(A)) £ 
K/[K*,K*], or Vet(A) = if A is not invertible. 

The Dieudonne determinant has a number of the desireable properties of the usual 
determinant, as proven in Dieudonne (1943): 

(1) Vet(AB) = Vet(A)Vet(B) for any A, B £ K" xn ; 

(2) Det(P) = 1 for any permutation matrix; 
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(3) Pet =Vet(A)T>et(B). 

\ob) 

Also note that if K has a degree function as above, then deg(Z>et(A)) is well defined, 
since all elements of the equivalence class of 7r(2?et(j4)) have the same degree (since the 
degree of all members of the commutator is zero. Gel'fand and Retakh (1991) show that 

Ssr(A) = |A| 11 |^ 11 )| 22 |A( 1 ^ 12 )| 33 |^ 123 ' 123 )|44 • • • \AV~ n -^~' n -V\ nn 
= \A\ n -Ser(A^), 

when all these quasideterminants are defined (or equivalently P is the identity in the 
Bruhat decomposition above), where ^(i- ■M - 'c) j s matrix A with rows 1 . . . k and 
columns 1 . . . k removed (keeping the original labelings of the remaining rows and columns) 
More generally, let R — (ri, . . . , r n ), C = (c%, . . . , c„) be permutations of {1, . . . , n}, 
let Rk = (ri, . . . , rfe), C/. = (ci, . . . , Cfc), and define A( Rk ' Gk > as the matrix A with rows 
n, . . . , Tk and columns c\, . . . ,Ck removed (where A( Ro ' Co ) = A). Define 

Sst r , c (A) = \A\ ruCl \A {RuCl) \r 2 , C2 \A {R ^ C2) \r^ 3 ■ \A^' c -% n , Cn 

= \A\ ruCl -8er RiC {A^^) (2.1) 
= \A\ ruC1 -\A^ c ^\ r2tC2 -Ser(A^ c ^). 

Fact 2.4 (Gelfand, Gclfand, Retakh, and Wilson (2005), Section 3.1). Let R,C be per- 
mutations of {!,... ,n\ and Rk,Ck defined as above. If \A^ Rk,Gk ' \ rk+1 ,c k+1 defined for 
k = . . . n — 1 , then 

Vet(A) = sign(i?) • sign(C) • w(5er RiC (A)). 

In other words, the Dieudonne determinant is essentially invariant of the order of the 
sequence of submatrices specified in (2.1). 

Theorem 2.5. Let A e K nxn be invertible, with degA^ < d. Then degPet(A) < nd. 



Proof. We proceed by induction on n. For n = 
predeterminants are 

8et 12 . 12 (A) 

SsT l2 ,2l(A) 

SeT 2 i,i 2 (A) 
Set 2 i,2i{A) 



1 this is clear. For n = 2, the possible 

A 12 A^ 2 A 21 )A 22 , 
A n A 2 iA 22 )A 2U 
A 22 A^A n )A 12 , 
A 21 A^A 12 )A 1U 



\A\nA 22 


= (A n 


\A\ 12 A 21 


= (A 12 


\A\ 21 A 22 


= (A 21 


\A\ 22 A n 


= (A 22 



at least one of which must be defined and non-zero, and all of which clearly have degree 
at most 2d. 

Now assume the theorem is true for matrices of dimension less than n. Choose r± , c\ £ 
{1, ...,n} such that |A| rijCl is non-zero and of minimal degree; that is deg^^^ < 
deg|A|fc^ for all k,£ such that \A\k,e defined and non-zero. The fact that |i| nA ^ 
implies that A^ Tl,Cl ^ is invertible, and we can continue this process recursively. Thus, 
let R — (ri, . . . , r n ) and C — (ci, . . . , c„) be permutations of {1, . . . , n} such that 
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|A' _R<,c ' i '| rj+1|Ci+1 ^ and deg | A^ Ri,Ci ^ |r i+ i,c i+ i is minimal over the degrees of non-zero, 
defined quasideterminants \A^ Ri ' Ci > \k.£, for < i < n. Now 

Ser R> c(A) = (|A| ri>cl • \A^ c \ 2tC2 ■ 5er(A^ c ^)) 

= Lo, -J^A rxk \A^\j k l A ec ^j ■ \A^% 2 , C2 -8er{A^ c ^) 

= A ruC1 -\A^ c % 2tC ^8er{A^ c ^) 

-53^*^.^17^ ■ \A^\ r2<C2 ■ 5er(A( R ^) 
k,e 

= A ruC1 ■ 5et{A^< c ^) 

-Y^Ar^^A^ ■ \A^\ r2<C2 • fer(A<** c »>), 
k,e 

where all sums are taken only over defined quasideterminants as above. Thus 

deg£>et(A) = Aeg8eT RyC (A) < max{d+ (n - l)d,2d+ (n- 2)d} < nd, 

using the induction hypothesis and the assumption that deg |A^ ri ' Cl ^ \ r2 ,c 2 1S chosen to be 
minimal. □ 

3. Degree bounds on matrices over F[d;a, S] 

Some well-known properties of f[d;a, 5] are worth recalling; see Ore (1933) for the 
original theory or Bronstein and Petkovsek (1994) for an algorithmic presentation. Given 
/, g G F[d;a, 8], there is a degree function (in d) which satisfies the usual properties: 
deg(fg) = deg / -I- deg g and deg(/ + g) < maxjdeg /, deg g}. We set deg = — oo. 

F[<9; cr, 8] is a left (and right) principal ideal ring, which implies the existence of a 
right (and left) division with remainder algorithm such that there exists unique q, r G 
F[d;a, 8] such that / = qg + r where deg(r) < deg (5). This allows for a right (and left) 
euclidean-like algorithm which shows the existence of a greatest common right divisor, 
h = gcrd(/, g), a polynomial of minimal degree (in d) such that / = uh and g = vh 
for it, v £ F[9;<7, 8]. The GCRD is unique up to a left multiple in F(f)\{0}, and there 
exist co-factors a,b G F[9;<7, 8} such that af + bg — gcrd(/, g). There also exists a least 
common left multiple lclm(/, g). Analogously there exists a greatest common left divisor, 
gcld(/, g), and least common right multiple, lcrm(/, g), both of which are unique up to 
a right multiple in F. From Ore (1933) we also have that 

deglclm(/, 3 ) =deg/ + deg.g-deggcrd(/,g), 
deg lcrm(/, g) = deg / + deg g - deg gcld(/, g) . 

It will be useful to work in the quotient skew field of F(d;a,8) of F[<?;<7, 8], and to 
extend the degree function deg appropriately. We first show that any element of F(<9; cr, 8) 
can be written as a standard fraction fg^ 1 , for f,g G F[d;a, 8] (and in particular, since 
F[<9;<7, 8] is non-commutative, we insist that g^ 1 is on the right). 
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Fact 3.1 (see Ore (1933), Section 3). Every element of F(d;o~,5) can be written as a 
standard fraction. 

The notion of degree extends naturally to F[d;a, S\ as follows. For /, g £ F[d;a,S], 
g ^ 0, wc define dcg(/g _1 ) = deg/ — degg for /, g e F[d;a,S]. The proof of the next 
lemma is left to the reader. 

Lemma 3.2. For f,g,u,v G F[d\cr, 8\, withg,v^=0, we have the following: 

(a) if fg^ 1 — to -1 then deg(/.g _1 ) = deg(uv _1 ); 

(b) deg((/ 9 - 1 ) • (uv- 1 )) = degifg-^ + degiuv-^J; 

( c ) deg(/5 _1 + uv~ l ) < max{deg(/g- 1 ),deg(uu" 1 )}; 

(d) deg((fg- 1 )- 1 ) = -deg(fg~ 1 ). 

In summary, the degree function on F(<9; tr, 5) meets the requirement of a degree func- 
tion on a skew field as in Section 2. 

3.1. Determinantal degree and unimodularity 

We now characterize a matrix being unimodular as those such that the degree of the 
Dicudonne determinant is zero. 

The following Fact forms the basis for the (theoretical) reduction of a matrix to Her- 
mite form. 

Fact 3.3 (Jacobson (1943), Section 3.7). Let a,b £ F[<9;er, 5], not both zero with g = 
gcrd(a, b), u, v € F[d;o~, 5] such that ua + vb = g, and s,( £ F[d\a 1 8] such that sa = 
—tb = lclm(a, b). Then 

W — i UV \ e F[d;a,S] 2x2 such that W 

\s t) \bj \n, 

and W us unimodular. 

Lemma 3.4. Let W be as in Fact 3.3. Then deg Vet W = 0. 

Proof. We may assume that gcrd(a, b) = g = 1, since the same matrix satisfies 
W(ag~ 1 , 6g _1 ) T = (1,0) T . Also assume both ab ^ (otherwise the lemma is trivial). 
Then 



ud| jaOj 1 1 v 
s t \ b 1 / \0t 



and Vet (W) ■ a = t mod [F[d; a, 8}*, F[d; a, 5]*}, 



so degdet W + deg a = degt. Since gcrd(a,6) = 1, from (3.1) we know deg a = degt, so 
degdet!F = 0. □ 

Theorem 3.5. LetU € F[d; a, 5} nxn . Then U is unimodular if and only if degVet U = 0. 

Proof. The Hermite form can be constructed (inefficiently) by continually zeroing col- 
umn entries with elementary unimodular matrices. Let w = (wi, . . . , w n ) T £ F[d; a, S] nx , 
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and i,j distinct indices such that Wi,Wj ^ 0. Let W be as in Fact 3.3 with (a, b) = 

(wi, Wj) T , and 

E = E(i, j; u, v, s, t) G F[<9; a, 5} nxn 
which is the identity matrix, except that 

En — 11, E{j — V, Eji — S, Ejj ^; 

where u, v, s, t G F[<9; a, 5] are as in Fact 3.3. Then E is unimodular, and Uw — (pi, . . . ,p n ) G 
F[<9; <t, <5]™ xl , with pi — gcrd(wi,Wj) and = 0. By Lemma 3.4 and the properties of 
Dicudonne determinants, AegDctE = 0. Finally, since the Hermite form of any unimod- 
ular matrix U G F[d; a, S] n n is the identity, any unimodular matrix is the product of a 
finite number of elementary unimodular matrices, and degPet {7 = 0. 

The converse is straightforward: if degl?et U — 0, then the Hermite form is the identity 
matrix, and the multiplier for the Hermite gives the inverse of U. □ 



3.2. Degree bounds on the Hermite form 

In this section we establish degree bounds on Hermite forms of matrices over F[<9; a, 8] 
and their unimodular transformation matrices. 

Theorem 3.6. Let A G F[<9; a, (5] nx ™ have full left row rank and entries of degree at most 
d and Hermite form H G F[<9; a, 5] n 71 . Then 

(a) The sum of the degrees of the diagonal entries of H has degree at most nd; 

(b) The sum of the degrees of the entries in any row of H has degree at most nd. 

Proof. Let V G F[<9; a, S] nxn be unimodular such that A = VH, whence T>et(A) = 
Vet(V)Vet(H). Therefore (a) follows from 

dcgVet(A) = dcgX>ct(ff) = degH u < nd. 

l<i<n 

Point (b) follows from the fact that each entry above the diagonal in the Hermite form 
has, by definition, degree smaller than the degree of the diagonal entry below it. □ 



We now show that all entries in H 1 have non-positive degrees. 

Lemma 3.7. Let H G F[<9; <r, S] nxn be of full left row rank and in Hermite form, and let 
J = H' 1 . Then deg J tj < for 1 < i,j < n. 



Proof. We consider the equation JH = 7, and note that J, like H is upper triangular. 
For each r G {1, . . . , n} we show by induction on c (for r < c < n) that deg J rc < 0. 

For the base case c — r, J rr H rr — 1, so deg J rr — — degH rr < 0. 

Assume now that r < c and deg J r g < for r < I < c. We need to show that 
deg J rc < 0. We know that 

^ JrtHf c = ^2 JrtHlc = 0. 
l<i<n r<£<c 

Since deg J r £ < for r < I < c and deg H cc > deg Hg c for r < £ < c, it must be the case 
that deg J rc < as well. □ 
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Theorem 3.8. Let A £ F[d; a, S] nxn be a invertible (over F(d;a,5)), whose entries all 
have degree at most d. Suppose A has Hermite form H £ F[<9; er, S] nxn , with UA = H 
and UV = I for U,V £ F[d;a,S] nxn . Then degV < d and deg[/< (n-l)degA 

Proof. Note that V = AH -1 , and by Lemma 3.7 all entries in H~ Y have non-positive 
degree. Thus d&gV < deg A. By Corollary 2.3, degU < (n — 1) deg A. □ 

4. Computing Hermite forms by linear systems over F[d;cr, S] 

In this section we present our polynomial-time algorithm to compute the Hermite 
form of a matrix over F[d;a, 6]. This generally follows the "linear systems" approach of 
Kaltofcn et al. (1987), and more specifically the refinements in Storjohann (1994) (for 
matrices over k[x] for a field k). We will need the tools for F[d;a, 5] we have developed 
in the previous sections. The method only directly constructs the matrix U such that 
H = U A. The Hermite form H can be found by performing the multiplication UA. 

Assume that Aij = J2o<k<d Aijkd k for Mjk G F. Let row(A,i) £ F[d;a,S] lxn be the 
iih row of A and define 

C(A) = l ]T b i -row(A,i):b 1 ,...,b n eF[d;a,S] 

[ l<i<n 

the left module of the row space of A. The following lemma is shown analogously to 
(Storjohann, 1994, §4.3.1, Lemma 4). 

Lemma 4.1. Let A G F[d;a,S] nxn be nonsingular, with Hermite form H. Let hi = 
deggHa for 1 < i < n. For v — (0, . . . , 0, ve, . . . , v n ) £ F[<9; er, S] xn , with deg ve < he, 
then if v £ C(A) we have vg = 0, and if vi ^ then v ^ £{A). 

The following Theorem is analogous to Storjohann (1994), §4.3.2, Lemma 5, with a 
different, slightly weaker degree bound. 

Theorem 4.2. Let A £ F[<9; er, S] nxn have full left row rank, with deggAij < d for 
1 < h J < n - Let (d\, . . . , d n ) be a given vector of non-negative integers. Let T be an nx n 
matrix with Tij = J2o<k<g^ijk9 k for unknowns tijk, where g > (n— l)<i-|-max i {(i i — hi}. 
Consider the system of eguations in tijk with constraints: 

(TA)i, iA = 1, for 1 < * < n, 
(TA) i>iik = 0, for k > d % , 
(TA)ij_k = 0, for i =/= j and k > dj 

By a solution for T we mean an assignment of variables tijk &ijk G F for some 
1 < z , j < n and < k < g such (4.1) holds. 

Let hi, . . . , h n £ F[d; a, 5] be the degrees of the diagonal entries of the Hermite form 
of A. The following statements about the above system hold: 

(i) If di > hi for 1 < i < n then there exist a solution for T; 



— diagonal entries are monic; 

— diagonal entry in row i has degree di\ 

— off diagonal entries have lower degree 
than diagonal entry in that column. 
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(ii) If there exists a positive integer £ < n such that di = hi for 1 < i < £ and dg < hg 

then there is no solution for T; 
(Hi) If di = hi for 1 < i < n then there is a unique solution for T such that G = TA is 
equal to the Hermite form of A under that solution. 

Proof. Let H e F[<9; a, S] nxn be the Hermite form of A and let U £ F[d; a, S] nxn be the 
unique unimodular matrix such that U A = H . 

To show (i), let D = di&g(d dl - h \ . . . , d d "- h ") e F[d; a, 5} nxn , and consider the 
equality DUA = DH. Let H* e F[d; a, 5} nxn be the Hermite form of DH and U* € 
F[<9; a, S] nxn the unimodular matrix such that U* DUA = H* . We construct H* from DH 
simply by reducing the entries above the diagonal (since it is already upper triangular). 
Thus U* is upper triangular, with ones on the diagonal, and deg U*j < (di — hi) — (dj — hj) 
for i < j. We claim T = U*DU is a solution to (4.1) . First note that the particular 
choice of D, together with the definition of H* ensure that the constraints of (4.1) are 
met. Furthermore, entries in the ith row of U* D have degree at most di — hi. By Theorem 
3.8, deg U < (n — l)d, whence degT < (n — l)d + max^c^ — hi} < q. 

To prove (ii), suppose by contradiction that there exists a nonnegative integer £ < n 
and a solution T such that degQ((TA)a) = di for 1 < i < I and deg 9 ((TA)«) < hg. Note 
that row(TA, £) = {{TA) e ,i, . . . , (TA)g, n ) is in C{A). First, if I = 1, then deg(TA)g A < h lt 
so by Lemma 4.1, (TA)i^ — 0, which is impossible since (4.1) ensures this entry is monic. 
Now assume I > 1. Then deg(TA)^i < hi (to satisfy (4.1)), and hence by Lemma 4.1, 
so (TA)g t i = 0. A simple induction shows that (TA)gj = for 1 < j < £. Now consider 
(TA)g^g, which has degree dg < hg by our assumption. Again by Lemma 4.1 (TA)g,g = 0, 
which (4.1) ensures is monic, a contradiction. 

If the conditions of (iii) hold, then by (i) there exists at least on solution for T. We 
can use an inductive proof similar to that used in our proof of (ii) to show that elements 
below the diagonal in TA are zero (i.e., that (TA)ij = for i > j). By the uniqueness of 
the Hermite form we must have TA = H. □ 



This theorem allows us to work with a partial order on the degree sequences. For 
any (hi, . . . , h n ), (di, . . . , d n ) G Z", we say that (hi, . . . , h n ) -< (di, . . . , d n ) if and only 
if hi < di for all 1 < i < n (and similarly define ^ for strict precedence). Thus, (4.1) 
has a solution if and only if (hi, . . . , h n ) -< (di, . . . , d n ) and this is unique if and only if 
(hi,...,h n ) = (di,...,h n ). 

We now embed the system (4.1) into a system of linear equations over F, with no Ore 
component. We embed F[d;a, S] into vectors over F via Tg : F[d;a,S] -4- F l+1 , with 

Tg(u + u x d + u 2 d 2 H h ugd^ 1 ) = (u , . . . ,u t ) e F l+1 . 

For g & F[d; a, 5] of degree d, u € F[d; a, S] of degree at most to, and assuming £> m + d, 
the equation ug — f can be realized by a matrix equation over F: 



/ n(g) \ 



(u , . . .,u m ) 




— (fo, ■■-,fe) 



.(«) vL(g) = n(f)- 



\Tt(d m g)J 
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Fixing di, . . . , d n € N as in Theorem 4.2, and setting g > (n — l)d + m&Xi{di — hi}, 
we can then study (4.1), as realized as (a subset of) the linear equations in the matrix 
equation over F: 

^(T n ) ••• T e (T ln )\ l^+ d (A u ) ■■■ ^+ d (A ln )\ /r e+d (G n ) ■■■ r e+d (G ln )\ 
\r e (T nl ) ••• T e {T nn )j \^+ d (A nl ) ••• nl+ d (A nn )J \r B+d {G nl ) ••• T e+d (G nn )j 



TGF" x <e+i)" 



AeF™(e+i)x(e+£i+i)r, 



GeF" x <e+d+i)« 



(4.2) 

We can ignore those equations whose component from G is unknown (i.e., not determined 
to be 1 or as in Theorem 4.2) because each such equation introduces no constraint on 
the solution. Thus, by Theorem 4.2, if we know d\, . . . , d n , this system will have a unique 
solution, from which we completely determine T . 

Example 4.3. Consider the following matrix in Q(z)[9;'] (the differential polynomials 
over Q(z)): 

' (z + 1) + d z + zd 

(z 2 + z) + zd z + l 2d 

^(— z — z 2 ) — zd zd zd J 

Assume for this example that we know the degrees of the entries in the Hcrmite form are 
(di,d2,dz) = (1, 0, 2). Then n = 3, and we can set g = 2, and have 



€Q(z)[d 
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We can remove columns from A for which there is an indeterminate in the corresponding 
column in G; these are coefficients of non-maximal degree in the column and are not 
specified in (4.1). Thus, we obtain the reduced system of equations 
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and can solve for 

2 z+1 

z 

2 z+1 

. 2z 2 +3z+2 
\ (z 2 +z+2)(2z+l) 

which corresponds to 

/ z+1 
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H = TA = 



{z + 1) + 8 
1 
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2z+l ' 



Z"+Z + 2; 



2z+l 



n 2z J +3z 2 -2z-5 a , a2 
U (z 2 +z+2)(2z+l)" + a , 



. /]3x3 



in Hermite form. 



We can now state our algorithm for computing the Hermite form give the degrees of 
the diagonal elements. 
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Algorithm Hermit eFormGivenDegrees 

Input: A e F[d;cr,6] nxn of full row rank, with (unknown) Hermite form H with diagonal 

degrees {hi,..., h n ) e N"; 
Input: (rfi, . . . , d n ) e N™, the posited degrees of diagonal entries of H 
Output: H e F[d; a, S] nxn if (di, . . . , d n ) — (hi, . . . , /i„), or a message that (di, . . . , d n ) 
is lexicographically smaller or larger than (hi, . . . , h n ); 
1: Let g = (n — l)d + max^ <ij 
2: Form the matrix equation = G as in (4.2) 

3: Remove all columns from G containing an indeterminate, and corresponding columns 
from A, to form the "reduced" linear system TA = G, where A and G are now 
matrices over F 

4: if rankvl < (n + l)g then 

5: return "(fti, . . . , h n ) ^ (di, . . . , d n ) v // System is underconstrained 
6: if TA — G has no solution then 

7: return "(hi, . . . , h n ) ^ (di, . . . , d n ) v // System is inconsistent 
8: Solve the system TA = G for f 
9: Construct T e F[d; a, <5]™ x ™ from f 
10: return H = TA and U = T 



From Theorem 3.6 we know that each entry in the Hermite form of A e F[d; a, <5]™ x ", 
with dcg a A < d, has degree at most nd. If the diagonal entries of A have degrees 
(hi, ... , h n ), then we know that 

(0, . . . , 0) -< (hi, ... , h n ) -< (nd, nd, . . . , nd). 

Algorithm HermiteFormGivenDegrees detects whether our choice of degree sequence is 
equal to, larger than, or not larger than or equal to the actual one. Thus, a simple 
component-wise binary search allows us to find the actual degree sequence (hi,..., h n ). 
That is, start by finding for the hi by executing HermiteFormGivenDegrees with degree 
sequence (di, nd, . . . , nd) for different values of di. This will require 0(log(nd)) attempts. 
Then search for hi using degree sequence 0(hi,di,nd, . . . ,nd) for different values of 
c?2, etc. It will require at most 0(n\og(nd)) attempts to find the entire correct degree 
sequence (hi, . . . ,h n ). 

Lemma 4.4. Given A G F[d;a,S] nxn of full row rank, where each entry has degree 
(in d) less than d, we can compute the Hermite form H € F{d;o~,5] nxn of A, and U G 
F[d;a,5] nxn such that U A = H. The algorithm requires us to call HermiteFormGivenDegrees 
0(n\og(nd)) times, with input A and varying degree sequences. 

For a first, general analysis of the complexity we will assume that operations in F have 
unit cost (and hence no coefficient growth is accounted for). To perform the rank test 
in Step 4, the inconsistency test in Step 6, and the equation solution in Step 8, we can 
simply do an LU decomposition of A using gaussian elimination. A has size n(g+ 1) x m, 
where n(g + 1) < m < n(g + d + 1), i.e., 0(n 2 d) x 0(n 2 d). Gaussian elimination can be 
accomplished with 0(n 6 d 3 ) operations over any field F (even if F is non-commutative). 
Combining this with Lemma 4.4 we obtain the following. 
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Theorem 4.5. Let A £ F[<9; cr, J]™*" have full row rank with entries of degree (in d) less 
than d. We can compute the Hermite form H £ F[<9; cr, 5] n n of A, and U £ F[<9; a, S] n n 
such that UA = H. The algorithm requires 0(n 7 d 3 log(nd)) operations in F. 

We next analyze our algorithm for computing the Hermite form of a matrix A £ 
k(z)[d; a, S] nxn over the field F = k(z), where k is a held and z an indeterminate. Without 
loss of generality A £ k[z] [d; cr, S) nxn by clearing denominators (which is a left-unimodular 
operation), but note that the Hermite form may still be in k(z)[d; cr, S] (see Example 
4.3). We will also assume for convenience that a(z) £ k[z] and deg<5(z) < 1. Thus 
dz = a(z)d + S(z) £ k[z][d] and the degree in z and d remains the unchanged. A more 
general analysis could follow similarly. 

Multiplying two polynomials in k[z] of degree at most m can be accomplished with 
0(M(m)) operations in k: M(m) — m 2 using standard arithmetic or M (to) = to log to log log to 
using fast arithmetic (Cantor and Kaltofen, 1991). We similarly assume that two integers 
with £ bits can be multiplied with 0(M(l)) bit operations. 

Theorem 4.6. Let A £ k[z] [d; cr, S] nxn have full row rank with entries of degree at most 
d in d, and of degree at most e in z. Let H £ k(z)[d;a,d] nxn be the Hermite form of A 
and U £ k(z)[d; cr, 5} nxn such that UA = H. 

(a) deg z Hij £ 0(n 2 de) and deg z f/y £ 0{n 2 de) for 1 < i,j < n. 

(b) We can compute the Hermite form H £ k(z)[d;a,S] nxn of A, andU £ k(z)[d; a, S] nxn 
such that UA = H with a deterministic algorithm that requires 0(n 7 d 3 log(nd) • 
M(n 2 de)) or 0~(n 9 d 3 e)^ , operations in k. 

(c) Assume k has at least 4n 2 de elements. We can compute the Hermite form H £ 
k(z)[d; cr, S] nxn of A, and U £ k(z)[d; cr, S] nxn such that UA = H using an algorithm 
that requires an expected number 0(n 7 d 3 log(nd) + n 7 d 3 e) of operations in k (using 
standard polynomial arithmetic). This algorithm is probabilistic of the Las Vegas 
type (never returning an incorrect answer). 

Proof. To show (a), recall that the matrix A £ F[z] is of size 0(n 2 d) x 0(n 2 d). Using 
Hadamard's bound and Cramer's rule, the numerators and denominators in T thus have 
degree at most n 2 de in z. H = UA has the same degree bound in z. 

To prove (b), we note that n 2 de is also a bound on the degrees (in z) of any reasonable 
algorithm for computing T. For example, we could simply compute modulo an irreducible 
polynomial in k[z] of higher degree than twice Hadamard's bound, and recover the entries 
in k(z) by rational recovery. The stated total cost follows. 

To show (c), we note that the tests for rank deficiency in Step 4, and inconsistency 
in Step 6, can be done by considering the equation T A — G mod (z — a) for a randomly 
chosen a from a subset of k of size at least 4n 2 de. This follows because the largest invariant 
factor w £ k[z] of A has degree at most n 2 de by Hadamard's bound (see part (a)), and 
the rank modulo (z — a) changes only if a is a root of w. By the Schwartz-Zippel Lemma 
(Schwartz, 1980) this happens with probability at most 1/4 for each choice of a (and 
this probability of error can be made exponentially smaller by repeating with different 
random choices). Thus, these tests require only 0(n e d 3 ) operations in k to perform, 



t We employ soft-Oh notation: for functions a and ip we say a £ 0~(ip) if <x 6 O (<p log c ip) for some 
constant c > 0. 
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correctly with high probability. During the binary search for the degree sequence we only 
perform these cheaper tests, requiring a total of 0(n 7 d 3 log(nd)) operations in k before 
finding the correct degree sequence. 

Once we have found the correct degree sequence, we employ Dixon's (1982) algorithm 
to solve the linear system over k(z) (this is the fastest known algorithm using standard 
matrix arithmetic, and is very effective in practice). This lifts the solution to the system 
modulo [z — a) 1 for i = 1, . . . , 2n 2 de, where a is a non-root of the (unknown) largest 
invariate factor of A (i.e., is such that rank^4 = rank A mod (z — a)). Computing the 
solution modulo (z — a) 2n de is sufficient to recover the solution in k(z) using rational 
function reconstruction, since both the numerator and denominator have degree less than 
n 2 de by part (a); see (von zur Gathen and Gerhard, 2003), Section 5.7. A random choice of 
a from a subset of k of size 4n 2 de is sufficient to obtain a non-zero of the largest invariant 
factor (and hence not change the dimension of the solution space) with probability at least 
1/4 by the Schwartz-Zippel Lemma. In the first step of Dixon's algorithm, we compute 
the LU-decomposition of A mod (z — a) using 0(n 6 d 3 ) operations in k. We then lift the 
solution to TA = G mod (z — a) 1 for i = } . . . , 2n 2 de. Each lifting step requires 0(n 5 d 2 ) 
operations in k, yielding a total cost of 0(n 7 d 3 e). □ 



Finally, we consider coefficient growth in Q of Ore polynomial rings over Q(z). For the 
computation, once we have constructed the matrix A, we can bound the coefficient-sizes 
in T directly using Hadamard-type bounds. We can then employ a Chinese remainder 
scheme to find the Hcrmitc form using the above algorithm (or any other method, for 
that matter) . For example, we could simply choose a single prime p with twice as many 
bits as the largest numerator or denominator in the solution to (4.2) and then compute 
modulo that prime, in 7L v \z\\ the rational coefficients of H can be recovered by integer 
rational reconstruction from their images in Z p (von zur Gathen and Gerhard, 2003, 
§5-7). 

However, for the purposes of analysis, it is interesting to see how big these coeffi- 
cients can grow. We consider matrices in A £ Z[z][d;a,S] nxn without loss of generality. 
For convenience in this analysis (though not in complete generality), we assume that 
deg z 8(z) < 1 and a(z) <E Z[z], so dz ~ a(z)d + S(z) € Z[z]. 

For a polynomial a — oq + a\z + ■ • • + a m z m G let Halloo = max^ |a,|. For 

f = fo{z) + fi(z)d + --- + f r (z)d r £Z{z}[d;a,5],\et H/IU = max* H/^. Define WA^ = 
maxy H-Ajj'Hoo. In equation (4.2), the entries in A have size at most 

\\A\\W =ma^max{\\A ij \\ 00 ,\\dA ij \\ 00 ,...,\\d e A ij \\ 00 } e Z. (4.3) 

Theorem 4.7. Let A € Z[z)[d;a,S] nxn be of full row rank and such that deg d (A) = d, 

deg z (A) < e and \\A\\ l £ < (5. Then the Hermite form H e Q(z)[d; a, 5} nxn and U € 
Q(z)[d; cr, 5} nxn such that UA = H satisfy 

log ||i?||oo, log Halloo € 0(n 2 d(loge + log/3 + logn + logd)). 

Proof. Entries in A are polynomials in Z[z] of degree at most e and coefficient size 
at most /3. Every minor of A, and hence each entry in the solution T, is bounded by 
Hadamard's bound, which in this case is 



((l + e)/3(n 2 d)) 



2 ,^0(n 2 d) 
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(see Giesbrecht (1993) Theorem 1.5 for height bounds on polynomial products). □ 



By performing all computations modulo an appropriately large prime (as discussed 
above), we immediately get the following. 

Corollary 4.8. Let A <= Z[z] [d; a, S] nxn have full row rank with entries of degree at 
most d in d, of degree at most e in z, and ||^4||cxi < ft (where g — 0(n 2 d)). Let H G 
Q(z)[d; cr, 6} nxn be the Hermite form of A and U € Q(z)[d; a, 6} nxn such that UA = H. 

We can compute the Hermite form H G Q(z)[d; a, 5} nxn of A, andU € Q(z) [d; a, 6} nxn 
such that UA = H , using an algorithm that requires an expected number 0((n 7 d s log(nd)+ 
n 7 d 3 e) ■ M(n 2 d(\oge + log/? + log n + log d))), or 0~(n 9 d 4 e log (3) bit operations. This al- 
gorithm is probabilistic of the Las Vegas type (never returning an incorrect answer). 

The following corollary summarizes this growth explicitly over two common rings, the 
differential polynomials Q(z)[d; r ], and the shift polynomial Q(z)[d;S]. 

Corollary 4.9. Let A E Z[z][d;a,8] nxn be of full row rank and such that deg d (A) = d. 
deg z {A) < e, H e Q(z)[d; cr, 6} nxn the Hermite form of A, and U € Q{z)[d; a, 6} nxn such 
that UA = H. For both the differential polynomials Q(z)[d; f ] (wjere a(z) = z, 5(z) = 1) 
and the shift polynomials Q(t)[3;<5>] (where o~(z) = z + 1, 5(z) = 0), we have 

log HC^IU, log Halloo £ 0~(n 2 d{e + log||A|| 00 )) 



Proof. To show this for differential polynomials, we note that for a = J2o<i<d a i( 2 )® 1 
Z[z][d;f], ~ * ' 

**= E (■) E «.w w ^. 

0<j<e 0<i<d 

where ai(z)^' is the jth derivative of di(z). Since only the first e derivatives of any ai 
are non-zero 

H^olU <i e - llalU -el 

and hence log ||^4.||^£? G 0(log WA^ + elog(n 2 d)) for q = 0(n 2 d). The result follows by 
Theorem 4.7. 

To show this for shift polynomials we note that for a = J2o<i<d a i( z )^ G ^[d, S], 

d e a= Yj a t (z + £)d\ 

0<i<d 

SO 

H^oHoo < ||a||oo2 e/2 r, 

and hence log ||^4||oo € 0(log ||^4.||oo + elog(n 2 d)) for q — n 2 d. Again, the result follows 
from Theorem 4.7. □ 
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